cty.info <- get_acs(geography = "county",
variables = c(pop="B01003_001",
medincome = "B19013_001"),
survey="acs5",
state = c("WA", "OR", "ID", "MT", "WY"),
geometry = TRUE, key = censkey, progress_bar=FALSE) %>%
select(., -moe) %>%
pivot_wider(
names_from = "variable",
values_from = "estimate"
)
p <- ggplot(data=cty.info) +
geom_sf(mapping=aes(fill=medincome))